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SUMMARY 


This manual describes the input data required to 
execute the computer program AIC / INT, a nemonic 
derived from the problem title "DEVELOPMENT AND 
APPLICATION OP SUPERSONIC UNSTEADY CONSISTENT AERO- 
DYNAMICS FOR INTERFERING PARALLEL WINGS". (Bell 
Aerospace Report No. 2471-941001). A finite element 
analysis is employed. 

The work was accomplished by Bell Aerospace 
Company under Contract NAS1-10880 with the National 
Aeronautics and Space Administration, Langley Research 
Center, Hampton, Virginia. 

A sample case is included that shows both input 
and output. 
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INTRODUCTION 


This manual describes the input data required to 
execute the computer program AIC / INT (aerodynamic 
influence coefficients with interference). Applicability 
is aided by the theoretical overview (see Ref. 1) and a 
sample case, including output. 

The purpose of the computer program is to generate 
aerodynamic forces for a pair of plane and interfering 
nearly parallel, non-coplanar wings at supersonic Mach 
numbers. A finite element technique has been employed. 
Planforms are described by triangular elements and dia- 
phragm regions are generated automatically. 

There are limitations to the program, in physical 
size (number of elements, grid points, etc.) and in 
analysis. The former is covered in this manual, the 
latter is found in Ref. 1. Of particular importance are 
the following 

a. Wake calculations are not included 

(e.g. in a wing-horizontal tail combination 
the influence of. the wake of the wing on the 
angle of attack of the horizontal stabilizer 
is ignored.) 

b. Results for an isolated wing are generated 
when placement of the wings is such that no 
interference occurs. (i.e. place a one element 
wing directly above the main wing at large 
separation. ) 


- 1 - 


SECTION 1 


FORMULATION 


The integral equations used for the computer program are 
presented here. A complete derivation of these equations will be 
found in Reference 1. 


1. Velocity Potential - Interferring Case (Subprogram TRINTX) 


The velocity potential influenced coefficient matrix is given by: 
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The initial position of each element is assumed to be the 
average of its Z-coordinates . Integrands in Eq. (1.1) are 

analytic in the fore Mach cone region. A Gaussian type 

/ 

numerical integration with 5 pivotal points is employed in 
all cases . 


2. Downwash Integrals (Subprogram DWASHY) 

The downwash integral, after performing an integration by 
parts and simplifying, becomes (See Reference 1) 
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( 1 . 2 ) 



where 



(1.3) 


(l.H) 


Since the singularities have been eliminated in these integrals, 
the limits of integration for the elements can extend up to the 
Mach cone. A typical element cut by the Mach cone is shown in 
Figure 1.1. Only partial integration is required here.. 
Numerical integration of the form 

r x T /l * r 

\ J Vi }(<?,*■) A A 

is to be performed. 

At a typical point, ^the upper limit Y u is. redefined as Y^ and 
the integration is stopped at the Mach cone. This technique is 
applied to all partial elements. 
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The input downwash terms are given by: 

1 kJ , > i » ii ■ 

w Pt . ° <j * 

where is the displacement and (<V<3* ) 

the slope in the stream direction of each node point. 


dj . Condensation of Velocity Potentials (Subprogram ASSPOT) 



where subscripts W and D denote the wing and diaphragm 
nodes. 
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6. Generalized Forces (Subprogram OUTPUT) 

From the definitions in Reference 1, 

( 1 . 12 ) 

( 1 . 13 ) 


- the d th deformation mode 

P = the nodal force distribution due to the deformation 
m 

mode X 
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Table 1.1 -Symbol Definitions 
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SECTION 2 


INPUT OVERVIEW 


The input package for this program expects to find required 
data arranged into specific groups. Each group of data, which 
may be one or more cards, is preceded by a LABEL card carrying 
an applicable name. These LABEL cards are used by the program 
as (1) an identifier for branching, and (2) a check that the 
data has been input. The input package supplied with this program 
will accept only the following LABEL data cards: 

* . 1. RUN 

2. TITLE 
3- SYSTEM 
H. LEDGE 
5- ELEM 
6. COORD 
7- MODE 

8 . END 

9. ENDATA 

10. FREQ 

LABEL data sections are subject to the following rules: 

Sections 1 and 2 are optional; 

Sections 4, 5, 6, and 7 may appear in any order; 

Section 3 must precede section 4, 5, 6, 7, 8; 

Section 8 must be the last data card to describe a wing; 
Section 9 must follow the second END card; 

Section 10 must follow ENDATA. 

Failure to observe these rules will terminate execution. 

Figure 2.1 gives- an overview of physical arrangement of the 
total data deck. 
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Figure 2.1 Input Overview 

Physical arrangement of data deck. Only the LABEL 
data header cards are shown in their relative locations. 
Note the absolute position of the FREQ card (s'). 
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Specific data formats correspond to each LABEL section. 

The formats and the LABEL section using same are presented here. 


TYPE FORMAT (---TYPE-.--) 


LABEL SECTION USING 


1 

2 

3 

4. 

5 

6 


20A4 

A4 , A2 j 14 
5X, 1215 

A4, 2X, 14, 1115 
5X, 15, 4E15.6 
10X, 4E15-6 


TITLE 

Label Cards Only 
SYSTEM (cardl); LEDGE 
ELEM (includes REPEAT) 
COORD 


MODE, FREQ 



SECTION 3 
INPUT FORMAT 


Input data is necessary for each subsonic/supersonic wing 
Therefore, to fully describe the problem, two complete (and 
independent) sets of data are required. The following pages 
describe each LABEL data section. in detail. 

Data is discussed and prepared under the following 
assumptions : 

1. Label data applies to each planform, one at 
a time (differences noted in the text). 

2. Wing 1 is chosen as the leading wing. Wing 2 

is normally above Wing 1 (i.e. Z 2 > 0) but this is 
not a requirement. 
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CARD COLUMN 


12 3 

! R • U N FORMAT: TYPE 2 

This card is optional. It is used to control partial or 
total execution of the entire program. This card may appear any- 
where before the end card. 

Case 1. Card Missing 

Program will process all data presented. Labeled 
sections will be accounted for as encountered. Problem 
size will be computed on the basis of system input data. 
All input data will be displayed in card image form in 
the order received. 

Execution will terminate after generation of the 
diaphragm grid point coordinates and element data. 

Case 2. Card Present 

All of Case 1 procedures plus total execution of 
the program 


CARD COLUMN 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

l T 

I 

T 

L 

E 





N 


FORMAT: TYPE 2 


This card is optional. 

If present, it must be followed by exactly N cards. Type 1 
format, containing title information. The information will be 
displayed exactly as presented, a useful aid in keeping track 
of output data. Maximum N = 20. 


This section may appear anywhere before the END card. 




SYSTEM 


Control Variables Definitions 


ITEM 1: 

ITEM 2: 

ITEM 3= 
ITEM 4: 

ITEM 5: 
ITEM 6: 


ITEM 7: 


ITEM 8: 


Total number of grid points on the wing. This 
number will be used to input COORD and MODE data. 

Total number of elements on the wing. This number 
will be used to input ELEM data. 

Total number of degrees of freedom for this system. 

Total number of grid points along the leading edge 
of the first sub-sonic crank of the wing. 

(See LEDGE section.) 

Number of equal divisions in the forward mach cone 
in the diaphragm region. (See Figure 3.1) 

Symmetry factor. Since only the right half of the 
wing is presented as physical data, the program 
must be told if the motion of the other half is 
symmetric (+i ) or antisymmetric (-1) 

Extrapolation of leading edge s wing 2 only. For wing 
1 leave blank. In situations 'shown in Figure 3 • 2 , 
the stagger of the second wing is such that the lead 
point on the upper wing lies within the forward and 
rearward mach cones of the lower wing. In this case, 
the leading edge is extrapolated from the lead point 
of wing 2 up to the intersection with the rearward mach 
cone of wing 1. The user must know if such extrapolation 
is required. 

Wing 1 leave blank. Wing 2 must supply a value 

1) of new points to be entered on the extrapolated 
portion of the leading edge. If Item 7 is zero or 
blank, Item 8 may also be neglected. 
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ITEM 9. Mach number, to be used for both wings 
(wing 2 value ignored) 

ITEM 10. Reference length. All grid point coordinates 
will be normalized with the wing 1 value. 

ITEM 11. Wing 1: This value used for truncation of a 

partial element near the mach cone 
during integration. A suggested value 
is 0.5 percent of the root chord. 

Wing 2: This value is used as the stagger distance, 

measured from lead point of wing 1 to 
lead point of wing 2. (the distance MP 
in Figure 3.2) The stagger distance is 
normalized by the reference length before 
it is used to position the second wing. 
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CARD COLUMN 



FORMAT -.TYPE 2 


This card is mandatory. 

It may be preceded only by optional cards. The data presented 
here pertains to the general properties of the system. It must be 
followed by 2 cards. 


Card 1 . Format : Type 3 

ITEM CARD COLUMN 


1 

NBR. of grid points on the wing 

6 

10 


planforra (NP) 






1 



11 

1? 

2 

NBR. of elements on the wing 








16 

20 

3 

NBR. of degrees of freedom 








21 

25 

4 

NBR. of leading edge grid points — 





1 



26 

30 

5 

NBR. of divisions in diaphragm — 








31 

35 

6 

Symmetry 





1 



36 . 

40 

7 

Extrapolation — 






8 

NBR. of extra grid points on dhord 

41 

45 


for extrapolation — — -- 







- 19 - 










Card 2. 


(System) Format: Type 6 









CARD COLUMN 


1 

2 

3 

4 

5 

UJ 

E 

D 

G 

1 E 1 


FORMAT: TYPE 2 


This card is mandatory. It may appear anywhere after SYSTEM 
section. It must always be present, even for super-sonic 
leading edges. 

The data presented here concerns the leading edge grid point 
numbers. Start with foremost point on the wing edge that is 
subsonic. Use 1 or more cards, type 3 format, to describe this 
data. The example below should clarify the input required. 



Case A. Both cranks sub-sonic 


SYSTEM, card 1, item 4=5 


LEDG] 

7 

—i 


LXJ 

2 

3 

4 

5 



Case B. Crank 1 is supersonic, crank 2 is subsonic. 

(Since the Mach line drawn at node 1 makes Crank 1 
supersonic, it has been shifted to node 3 for clarity,) 
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CARD COLUMN 


1 

2 

3 

4 

5 

E 

0 

0 

1 R 

d| 


FORMAT: TYPE 2 


This card is mandatory. 

It may appear anywhere after SYSTEM. There must follow exactly NP 
(SYSTEM, card 1, item 1) cards, type 5 FORMAT, containing grid point 
coordinate data for the planform in questions. Grid point numbers 
must be generated consecutively, beginning with the number 1 (one). 

NOTE: Remember that this data will be divided by reference 

length supplied on card 2, system.. 


Grid point nbr 


X-coordinate - 


Y-coordinate - 


Z-coordinate - 


6 7 8 9 10 


11 25 



- 22 - 








CARD COLUMN 


1 

2 

3 

4 

1 E 

L 

E 

1 M 1 


FORMAT: TYPE 2 


This card is mandatory. 


It may appear anywhere after SYSTEM. The total number of cards to 
follow must be equal to SYSTEM, card 1, item 3* The cards are 
type 3 format. 


Each card contains 4 numbers as follows: 


Element Number 

8 

9 

10 





13 

14 

15 

Grid point # 1 



1 


18 

19 

20 

Grid point #2 





23 

24 

25 

Grid point #3 





NOTE: 

Element numbers must be generated consecutively, beginning 
With the number 1 (one). Therefore, the total number 
of elements on the planform (SYSTEM, card 1, item 2) must 
be equal to the highest element number. 
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An alternate method of entering element grid point 
data is available to the user in the form of a REPEAT option. 
This option can reduce considerably the amount of input data 
required to define the structure. (Refer to FIGURE 3>3 and 
Table 3 . 1 ) 

This technique will be useful only when the grid points 
of the next higher (numerically) elements (say 1 to N) have 
their grid point coordinates increased by a common factor. 

For this case, one need only define a first element and supply 
a repeat "card (type 4 FORMAT) as follows: 

CARD COLUMN 


Repeat how many times? 



Increase grid point #1 by 


Increase grid point # 2 by 



Increase grid point #3 by 23 24 25 
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Figure 3.3 Sample Planform of 15 elements 
and 14 gridpoints. Element 
numbers are circled. 
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Pull , Input 

ELEM 


1 

6 

2 

1 

2 

7 

3 

2 

3 

8 

4 

3 


9 

5 

4 

5 

7 

6 

2 

6 

8 

7 

3 

7 

9 

8 

4 

8 

10 

7 

6 

9 

11 

8 

7 

10 

12 

9 

8 

12 

12 

11 

8 

13 

13 

11 

10 

14 

14 

12 

11 

15 ' 

14 

13 

11 

11 

11 

10 

7 


Mixed Input 

ELEM 


1 

6 

2 

1 

REPEAT . 3 

1 

1 

1 

5 

7 

6 

2 

REPEAT 2 

1 

1 

1 

8 • 

10 

7 

6 

REPEAT 2 

1 

1 

1 

11 

11 

10 

•-f 

1 

12 

12 

11 

8 

13 

13 

11 

10 

14 

14 

12 

11 

15 

14 

13 

11 


Table 3*1 ELEMent grid points Data for 
the Sample structure in 
Figure 3.3 
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12 3 4 


M 

0 

D 

E 


FORMAT: TYPE 2 


This card is mandatory. It may appear anywhere after system. 

Data presented here refers to the mode shapes applicable to this 
test/planform. The data is considered an array, 

card position implies grid point number, and column position the 
relative mode shape. Therefore, there must be NP (SYSTEM, card 1, 
item 1) cards of type 6 format. 


Note : 


Mode shape data is used as received. No normalization process 
occurs. Therefore, the user must be careful to insert the correct 
mode shape data. 

I . 
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CARD COLUMN 



FORMAT: TYPE 2 


This card is mandatory. 

It must physically appear as the last card of a wing (planform) 
description. 


CARD COLUMN 


1 

2 

3 

4 

5 

6 

E 

| N 

D 

A 

T 

A | 



FORMAT .-TYPE 2 


This card is mandatory. 

It must appear immediately after the second END card, i.e. at the 
end of the label data describing wing'2. If this. card is missing, 
execution is terminated before any data is processed. (I.E. The 
program strikes an END-OF-FILE condition, and reacts accordingly). 
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CARD COLUMN 



This card is mandatory. 


The numerical value ( K ) placed on this card is called 
the "reduced frequency" and is determined by 


where 

Oj 

= angular flow frequency 


L 

= reference length (SYSTEM 


XT 

= flight speed. 


Although H is normally considered to be complex, the real 
part is taken as zero and the user enters only the magnitude 
of the complex part as a floating point number in columns 
11 thru 25. To compute generalized forces for more than 
one frequency, simply enter additional cards. Each card 
contains only one frequency. No limit as to the number of 
cards allowed. 

Place all frequency cards immediately after "ENDATA" label 
card. 






SECTION 4 


SAMPLE PROBLEM/DATA/RESULTS 


The purpose of this section is to demonstrate graphically 
how the user describes a planform, the actual data required, and 
the results generated by the computer. Figure 4.1 shows a simple 
DELTA-WING planform consisting of 4-gridpoints and 2-elements. The 
broken line elements and gridpoints are generated interally, but 
the user should always have an idea of what this ’’DIAPHRAGM’’ region 
looks like before data preparation begins. The angle ©-(<.2,1,6) 
is computed thru SINE ©= 1 / (MACH NBR). Gridpoints are indicated 
by numbers. Elements are the circled numbers. Note that the 
coordinate system is a normal right-handed system if this page is 
turned 90-degrees. The CHORD of this planform lies on the X-axis, 
with gridpoint 1 at the origin. 

Should questions arise during the study of this example, 
re-read the LABEL . sect ion covering the question. 



^ 4 

Figure 4.1 Sample Planform Showing Basic Structure 
and Program Generated Elements 
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The printout from any run will have the same form as that 
presented here. (Note: this printout has been compressed to 

preserve continuity. The content remains unchanged.) Under conditions 
of supersonic wings, no diaphragm is generated and therefore no 
printout occurs for that section of the run. The following 
discussion will serve as a guide to understanding all results. 


ITEM 

DESCRIPTION 

PAGES 

A. 

An echo-check of the entire data deck 

33 

B. 

Process data, display after assignment to 
internal variables. Successful input pro- 
cessing yields a +++END CARD ENCOUNTERED+++ 

34-37 

C. 

Generation of diaphragm 

38-39 

D. 

Display input downwash 

40 

E. 

Formal output begins at TITLE 

41 . 

P. 

Wing 1, potentials, pressures, generalized 
forces 

42 

G. 

Wing 2, potentials, pressures, generalized 
forces 

43 

H. 

Repeat D thru G for each additional 
frequency 

— 

I. 

Job completed. Message END OF JOB 

44 


Abnormal termination of the program is defined by a) a 
full page message ERROR HALT, b) sufficient messages to indicate 
the reason, c) under certain conditions a core dump will be 
provided. Poor input data will normally be the problem. 
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The problem results (items P, G above are formally 
defined as follows: 




1. Velocity potentials (j) 

<p = - $ 


/*S> 

2. Pressure (in terms of nodal forces) p 

P = - 4 ib Ll P 


3. Generalized forces Q . . (complete wing) 


Qa - 


where 


-- M P v 


% ' t r 


is the dynamic pressure 


This completes the description of the computer printout 
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ELEMENT GRID POINTS 
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element grid points 
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DIAPHRAGM ELEMENTS WEPE 
ADDEO AS FOLLOWS. 









DIAPHRAGM ELE MEM $ kEPE 
ADDED AS FOLLOWS. 
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• • A I C / INT.. ..TEST CASE.. TEST0CC4 

THIS EXAPPLE IS USED TO OE HOMS TP ATE THE INPUT DATA TESTOOCD 

KCCUIKED TO EXECUTE THE *AIC / INT* PROGRAM. TEST0006 

BOTH WINGS ARE IDENTICAL, AND THE UPPER WING LIES TEST0007 
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COL 1 COL 2 

ROW REAL IMAG REAL { HAG 

*************************** *** ********* ***************************** 

1 0.113371E 00 0.2915636 CO 0.152803E 01 -0.3821196 00 

2 0.3808116-01 0.1063936 00 0. 555082E 00 -0.113576600 


PRtH Ndk.. 0.115JI 01 
FR6CU6NCY.. O.A000T 00 





COL 1 CCL 2 

ROW P.6AL; IPAG Pf AL I PAG 

1 0. 1 133746 00 0.2915636 00 0. 1528036 01 -0.3621A9E 00 

2 0.3808 146-01 0.1063936 00 0.5550836 00 -0.113576E 00 





SECTION 5 

RESTART / PUNCHED OUTPUT 

The AIC/INT program has an additional analytic feature 
called "RESTART". What the word restart really means is that 
the user may alter MODE shape data and generate new aerodynamic 
coefficients without making all of the initial calculations 
performed on the initial run. An assumption is made that the 
data from the initial run was permanently stored. 

Punched output (generalized forces and pressures) is also 
an option. Each block of data punched contains a descriptive lead 
card containing the type (forces or pressures), wing number and 
mach number. 

To exercise either of : the above options, code in columns • 

9 and 10 of the run card as follows: 

+1 = punch generalized forces 

+2 = punch generalized forces and pressures 

-1 = punch generalized forces and restart 

-2 = punch generalized forces, pressures and restart 

The format of the punched output is (4E15.6), 


-4 5- 


SECTION 6 


HINTS ON USING AIC/INT 


The following comments apply to every planform the user 
wishes to test. 

a. Create a scale drawing of the planform. 

b. Create finite element pattern by 

1. selecting gridpoints first and draw 
connecting lines) 

or 

2. draw lines first and let intersections 
be gridpoints . 

c. Label grid points in increasing order, from 
0 to +x, 0 to +y. 

d. Label elements as in c. above, numbering all 
elements in increasing order starting with 1 (one). 


e. Prepare ELEM data by choosing gridpoint #1 
to have the largest x value of the 3 points 
involved. (Any method is acceptable.) 

f. Prepare COORD data. 

g. Prepare remainder of data. 

h. Submit data to execution but leave out RUN card. 

i. Check the printout of the input data carefully for 
errors. Compare the grid point data with the scale 
drawing of item a. 


j. Prepare frequency cards, enter after ENDATA card, 
include a RUN card and begin execution. 
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SECTION 7 


COMPUTER PROGRAM FLOW 
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1 


Assembly 

CALL 


elem. 0 to 

ASSPOT 

— — — — 

master (fi and 



condense 


CALL 

ASSDWN 


Assembly 
downwash pots 
compute system 
analysis 




Display 

input downwash 

CALL 


DISPLA 


array 



Display 
final results. 
Wing 1 


Display 
final results 
Wing 2 


CALL 

ENDJOB 


STOP 


—s 

) 
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.-SECTION 8 
SIZE LIMITATIONS 

At delivery date this program contains the following 
limitations (applicable to planform) : 


1. 

Nbr . 

gridpoints on wing 

4 - 

Ho 

2 . 

Nbr . 

gridpoints on diaphragm 


15 

3* 

Nbr. 

degrees of freedom 


10 

4. 

Nbr . 

leading edge gridpoints 

/ 

20 

5 • 

Nbr . 

elements (wing + diaphragm) 

Z. 

100 

6. 

Nbr . 

title cards 


20 


These limitations are subject to change. 
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